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Our account of the problem of the classical limit of quantum mechanics involves two elements. 
The first one is self-induced decoherence, conceived as a process that depends on the own dynamics 
of a closed quantum system governed by a Hamiltonian with continuous spectrum; the study of 
' decoherence is addressed by means of a formalism used to give meaning to the van Hove states with 

diagonal singularities. The second element is macroscopicity represented by the limit ft — > 0: when 
the macroscopic limit is applied to the Wigner transformation of the diagonal state resulting from 
decoherence, the description of the quantum system becomes equivalent to the description of an 
ensemble of classical trajectories on phase space weighted by their corresponding probabilities. 
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I. INTRODUCTION 



> 

The problem of the classical limit of quantum mechanics has been a point of debate since the birth of the theory. 
Although this problem is usually addressed in the context of measurement, it can be analyzed from a more general 
' point of view, in terms of how the classical world arises from an underlying quantum reality, independently of whether 
| there is a measurement involved or not. Of course, the problem of the classical limit relies on the assumption that, 
if quantum mechanics is correct, then its results must reproduce the results of classical mechanics in the appropriate 
limit. 

In the old days of the theory, Heisenberg and Bohr among others conceived the classical limit of quantum mechanics 
by analogy with the classical limit of special relativity: h — > in quantum mechanics should play the same role as 
f3 — > in special relativity. This assumption was considered by Einstein as an oversimplification since, while relativity 
and classical mechanics have the same deterministic structure, quantum mechanics has a probabilistic structure. 
Nevertheless, since those days it has been usually claimed that classical mechanics can be recovered as a limiting case 
of quantum mechanics when H — > 0. This assumption led to correct results when the classical limit was conceived in 
the following way: 

X 



Qm| quantization 1 CM 

(_ classical limit = n — > J 

where QM and CM stand for quantum mechanics and classical mechanics respectively. In this schema, the first step is 
to quantize a classical system, e.g., by means of the Weyl transformation, in order to obtain the corresponding quantum 
system (at present, quantization is also called "deformation"). Then, the original classical system is recovered by 
applying the inverse Weyl transformation, i.e., the Wigner transformation 1 to the quantum system previously obtained 
and by taking the limit h — ► 0. It is quite clear that this method is completely circular to the extent that it only 
recovers the classical system originally proposed. 

When the theoretical structure of quantum mechanics finally lost its classical origin, the problem of the classical 
limit acquired a new formulation that became the traditional one: 



qm{" 



classical limit = h — > > CM (2) 



1 Historically, Weyl proposed his transformation as a quantization method. Later and independently, Wigner proposed a trans- 
formation that mapped quantum states into classical density functions. Finally, Moyal proved that the Wigner transformation 
was equivalent to the inverse Weyl transformation. 
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This means that, no matter how the original quantum system was described, a classical system should be obtained 
via the Wigner transformation when ft — > 0. However, this way of conceiving the problem of the classical limit leads, 
at least, to three problems: 

1. In general, the Wigner state function p{4>) (where 4> = {liP) is a point in phase space) is not non-negatively 
defined 2 ; as a result, it cannot be interpreted as a probability distribution. 

2. Only Hamiltonians of degree < 2 in p and q yield to Hamiltonian fluxes that maintain the deformation invariant 
(or covariant) (see [1], [2]). In fact, only in these cases p\pi{t) — p\(t)pi(t) after performing the Wigner 
transformation. 

3. In some cases, factors of the form ft -1 may appear in the Wigner state function due to the features of the Wigner 
transformation (see [1], [3]). In these cases, the limit ft — > of the Wigner function is singular. 

In this paper we will follow a well known trend in contemporary physics, according to which the classical limit 
must not be conceived as a mere consequence of a limiting procedure, but as a result of a physical process. From this 
perspective, the explanation of the emergence of the classical world from the underlying quantum realm involves two 
steps: the first one consists in explaining the physical phenomenon of dccohercnce, and the second one consists in 
taking the macroscopic limit ft — > 0. However, we will move away from the mainstream position with respect to the 
explanation of decoherence: the aim of this paper is to obtain the classical limit of quantum mechanics on the basis of 
the self-induced approach to decoherence, such as it was presented in paper [4] and discussed in depth in paper [5] . In 
contrast to the traditional einselection approach [6], from the self-induced perspective decoherence does not require 
the openness of the system and its interaction with the environment: a single closed system can decohere when it 
has continuous spectrum. We will show that, in this new scenario, the classical limit is described by the following 
diagram: 



„ J decoherence Boolean QM macro scopicity = h — > 1 ,~, 01 . , toX 

(jlvl < _ > > CSM (3) 

L classical limit J 

Self-induced decoherence transforms quantum mechanics into a Boolean quantum mechanics where the interference 
terms that preclude classicality have vanished. Macroscopicity, expressed by the limit ft — > 3 , turns Boolean quantum 
mechanics into classical statistical mechanics (CSM) in phase space. According to this view, the classical limit 
of quantum mechanics is not classical mechanics but classical statistical mechanics, and it requires two physical 
conditions: decoherence and macroscopicity. In other words, in order to behave classically a quantum system must 
have decohered and must be macroscopic enough: each one of these conditions alone is necessary but not sufficient 
for its classical behavior. Furthermore, we will show how and under what conditions this explanation overcomes the 
three problems that arise from the traditional way of conceiving the classical limit. 

This paper is organized as follows. In Section II we present the formalism for observables and states, necessary 
for developing our program. Section III is devoted to explain the self-induced approach to decoherence: decoherence 
in energy and in the remaining variables arc considered. In Section IV we study the operation known as Wigner 
transformation and its application to observables and states. In Section V we show how the classical limit leads to 
classical statistical mechanics when the macroscopic limit is applied to the Wigner transformation of the quantum 
state resulting from decoherence. In Section VI we discuss the physical meaning of the results just obtained, arguing 
that classicality must be understood as an emergent property that objectively arises from an underlying quantum 
mechanical realm. Finally, in Section VII we draw our conclusions. Two appendices complete the paper. 



II. FORMALISM FOR OBSERVABLES AND STATES 



The formalism for observables and states used in the present paper is inspired by the formalism introduced by 
Antoniou et al. ( [7], [8]) which, in turn, is based on the works of van Hove [9]. In several papers ( [10], [11], [12], 



2 A function p(<j>) is non-negative if and only if p(0) > 0. 

3 It is quite clear that it is not possible to set the value of ft equal to 0, since it is not a dimensionless parameter but an 
universal constant. This means that, strictly speaking, the macroscopic limit is ft/5 — > 0, where S is the characteristic action 
of the system: this is a factual limit which represents realistic situations where S » ft. 
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[13], [14], [15]) we applied different versions of this formalism to the study of the properties of quantum systems 
with continuous spectrum. In particular, the formalism was used in paper [4] for explaining decohcrcncc. In order 
to simplify the notation, here we will study a simplified model where all the observables have continuous spectrum 
(cases where all observables except H have discrete spectrum will be considered in the footnotes): this will allow 
us to improve the mathematical basis of our approach without a proliferation of indices that would not introduce 
conceptual advantages. 



A. Quantum operator algebra 

Let us consider a system with a complete set of commuting observables (CSCO) {H, 0\, On} where H has 
a continuous spectrum < uj < oo and, for the sake of simplicity the Oi, i = 1,2,..., TV, have also continuous 
spectra 4 . We will assume that the observables H, 0\, On are Weyl observables, i.e., that they come from the Weyl 
transformation of classical observables. In order to simplify the notation we will use {H, O} to denote the CSCO 
{H, Oi, On}- The generalized eigenbasis of {H, O} is {\uj, o)}, where uj and o satisfy: 

H \uj, o) = iv \ui, o) and O \ui, 6) = o \ui, 6) (4) 

Then, H and O can be expressed as: 

poc p poc p 

H = / uj \u),o)(u),o\ duo do , O = / / o \uj, o)(uj, o\ duo do (5) 

JO Jo Jo Jo 

In addition to H and O, there are additional observables that may or may not commute with H and O. Then, a 
generic observable A has the following form: 

pOO pOO P P 

A= / / / / A(u, uj' , o, d) \u, o)(uj', o'\ duo did dodo (6) 
Jo Jo Jo Jo' 

where A(w, uj', o, o') could be, in principle, a distributional kernel. However, we will not work with the set of all the 
possible observables of the system, but only with a subset of it. The condition that defines this subset is given by the 
choice of the kernel A(u, uj' ' , o, o'), which it is usually taken to be ( [7], [8], [12]): 

A(u>, uj', o, o') = A(uj, o, o') 5(u> - uj') + A(uj, uj', o, o') (7) 

where A(u),o, o') and A(u),u)', o, o') arc sufficiently regular functions (see [13] for details). Then, we will work with 
observables whose generic form is: 

A= / / / A(uj, o, o') \uj, o, o') duo dodo' + / / / / A(u>,u>',o,o') \uj,uj' ,o,o') dioduj' dodo' (8) 

JO Jo Jo 1 JO JO Jo Jo' 

where we have introduced \uj,o, o') = \u),o)(u),o'\ and \u),u',o,o') = \ui, o)(ui', o'\. With the condition {uj, o\uj' , o') = 
5(uj — uj') S(o — o'), the set of the operators of the form (8) is an algebra A, and the observables are the self-adjoint 
elements of A [12], [13] 5 . 

The first term of eq.(8) represents the observables that commute with those of the CSCO {H, O}, and it will be 
called the singular component As of A; the second term of cq.(8) will be called the regular component Ar of A: 



(9) 



f P POO P P POO POO 

As -— / A(uj, o, o') \uj, o, o') duj dodo' ; Ar := / / / / A(uj, uj' , o, o') \uj, uj' , o, o') duj duj' dodo' 

Jo Jo' JO Jo Jo' JO JO 

The operators As and A R form the singular and regular algebras As and Ar respectively: A s S -4s and A R G Ar. 
The generalized eigenbases of As and Ar are {\uo, o, o')} and {\uo, uj' , o, o')} since they span the algebras As and Ar 
respectively. Note that these two algebras have a trivial intersection and that the algebra A is the direct sum of both: 
A = As © Ar [13]. 



4 The continuous spectrum is relevant for the classical limit since, in the limit h — > (precisely, the high quantum number 
limit), many discrete spectra become continuous. Spectra with continuous and discrete parts are studied in papers [4] and [15]. 

5 Although we will work with a subset of all the possible observables of the system, the physical generality of the self-induced 
approach to decoherence relies on the fact that the coordinates of the observables not belonging to A in the generalized 
eigenbasis of {H, O}, being singular, cannot be measured in laboratory and, therefore, they must always be approximated by 
their averaged counterparts (for a full argument, see [5]) 
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B. States as linear functionals 



States are continuous linear functionals on the algebra A defined as above. Let A* be the dual space of A (the 
vector space of all linear continuous functionals on A). In our notation, the action of the functional p G A* onto the 
operator A e A is denoted as With this notation we define [ui, o, o'\ and (ui,w', o, o'\ as: 

(u>,o, o'\A) = A{u>, o, o') , (u,u',o, o'\A) = A(u,u',o, o') (10) 

for all A € A. It can be shown that (u>, o, o'\ and (u>, uj' , o, o'\ are in A* for all values of uj, u>', o, o' [13]. In addition, if 
Ag is the dual of As and A* R is the dual of Ar, it can be shown that (u>, o, o'\ e .4J and {u>,ui' , o, o'\ S *4Jj [12], [13], 
and that the following relations hold: 

(uj, o, s\u>', o', s') = 5(ui — u)')5(o — o')S(s — s') , (ui, a, o, s\w' , a', d , s') = 8(u> — u>')5(cr — a')S(o — o')S(s — s') , 



(oj, a, \ui', a', o', s') — (u, a, o, s\u>' , a') = (11) 
An element of the dual A* can be expressed as: 

pOO p p pOO pOO p p 

p= / / p(w, o, o') (uj, o, o'| duj dodo 1 + / / / / p(lu,lu' ,o,o') {uj,lo' ,o,o'\ du duo' dodo' (12) 

JO Jo Jo' Jo Jo Jo Jo' 

Again, we will call the first term of eq.(12) the singular component ps of p, and the second term of eq.(12) the regular 
component pn of p, where ps £ A* s and pa G A* R : 

ps '■= / / p( w i o, o') (u, o, o'| a\xj dodo' ; pr .= / / / p(ui, uj' , o, o') (w, u) 1 , o, o'\ dui du' dodo' (13) 
Jo Jo Jo' Jo Jo Jo Jo' 

The action of the functional p on the operator A is given by: 

pOC p p pOO pOC p p 

(p\A) = / / / p{u), o, o') A(u), o, o') duo dodo' + / / / / p(lo,lo' , o, a!) A(lo,lo' , o, o') duo duo' dodo' (14) 

Jo Jo Jo' Jo Jo Jo Jo' 

It is interesting to remark that, although A(u, o, o') and A(u, u', o, o') must be regular, well behaved functions (poly- 
nomials and Schwartz functions, see [13]), this is not the case of p(u, o, o') and p(uj 7 uj', o, o'), which may be singular. 
For instance, the functional (w, s, s'\ can be written in the form (12) with p(u>, o, o') — 5(w — uj) 5(s — o) 5(s' — o') and 
p(u>, uj', o, o') = 0. 

The condition of positivity for a functional / means that, if / € A* and A 6 A, then f(A^A) > 0, where A^ is the 
adjoint of A. In our case we will require positivity to ps, and this implies that: 

p(uj,o,o') >0. (15) 

The condition of normalization for a functional / means that the identity I must be an element of the algebra and 
that /(/) = 1. In our case, we will normalize only ps 6 : 



pOC p pOO p p 

1=1 / \uj,o){uj,o\ dujdo (ps\I) = / / / p(u>,o,o') du; dodo' = 1 

JO Jo Jo Jo Jo' 



(16) 



Finally, note that „4* = A* s © A* R [12], and that {(w,o, o'\} is a generalized basis of Ag and {(u),u)', o, o'\} is a 
generalized basis for A* R . 



We require positivity and normalization for ps since, as we will see, it is the only component of p that remains after 
decoherence. In addition, (pn\I) = for any pn 6 Ar. 
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III. SELF-INDUCED DECOHERENCE 



A. Decoherence in energy 

Let us now consider the time evolution of the system. Since p is a functional, its time evolution in the Schrodingcr 
picture cannot be directly computed by means of the Liouville-von Neumann equation. Nevertheless, this equation 
also describes the time evolution of the observable A in the Heisenberg picture: 

ihd t A = -[H, A] = hA => A{t) = exp(itL/ft) A(0) (17) 

where H is the Hamiltonian that governs the time evolution, and L is the Liouville operator associated to the 
Hamiltonian H: hA = — [H, A]. Once the time evolution of A has been computed, the time evolution of p can be 
obtained by means of the duality formula: 

(p\ exp(-ith/h) A) = (cxp(ith/h)p\A) (18) 

This equation gives the time evolution of p, which satisfies the Liouville-von Neumann equation: 

ihdtp = [H, p] = -hp => p(t) = cxp(i<L/%(0) (19) 

In order to follow this strategy in our case, we begin by applying the Liouville-von Neumann evolution equation 
to the generalized basis {\ui,o,o') = \u,o)(u),o'\, \uj,uj' ,o,o') = \ui, o)(ui', o'|}. Since H\uj,o,o') = uj \uj,o,o'), we have 
that: 

h\u>,o,o') = -H\u),o)(u>,o'\ + \uj,o)(uj,o'\H = -(u> - u>)\u>,o,o') = (20) 
L|w,w', o,o') = -H\uj,o)(uj',o'\ + \uj,o)(uj',o'\H = -(uj-uj')\uj,uj',o,o') (21) 

This means that the generalized basis {\u>, o, o'), \u, w' , o, o')} is an eigenbasis of the operator L. Moreover, h\u), o, o') = 
implies that not only the \u, o, o'), but also all the singular operators A$ G As are time invariant (e~ ltL / h As = As), 
since: 



p OO p p pOO p p 

hA s =h / / A(u),o,o') \ui,o,o')dLjdodo' = / / / A(w, o, o') (L \u>, o, o')) dui dodo' = 

^0 Jo Jo' JO Jo Jo' 

Therefore, for any A e A, 1LA = 1LA S + ~LA R = hA R . Moreover, due to eq.(21) ~LA R e Ar. 
From eq.(21), it can be obtained: 



(22) 



e~ uh/h \lo, uj', o, o') = e «(---')M | W) 0? ') ( 23 ) 
and, hence, for any Ar £ Ar we have: 



p p p OO p OO 

e ^tL/h AR = / / / / A(uj,uj' ',o,o')e itl - UJ - UJ '^ h \uj,uj l ',o,o')dujduj' dodo' (24) 

Jo Jo' Jo Jo 

Then, for any A £ A, we obtain the following time evolution in the Heisenberg picture: 

/■OO /■ r /"OO rOO /• r 

e - it h/h A= / / / A{Lu,o)\uj,o,o')dujdodo' + / / / / A(u,uj',o,o')e lt(ul - ul ">/ h \uj,uj',o,o')dujduj' dodo' 

Jo Jo Jo' Jo Jo Jo Jo' 

(25) 

A similar situation arises when we consider the time evolution of the states. Since \ui,o, o') is time invariant, the 
duality formula tells us that the functionals (uj, o, o'\ are time invariant; therefore, all the singular functionals ps G A* s 
are also time invariant. 

Now we can compute (p\A(t)) which, by the duality formula (18), is equal to (p(t)\A): 



( P \A(t)) = (p(t)\A) 



/ / / p(u), o, o') A(u), o, o') duo dodo' + 

Jo Jo Jo' 



/ / / / p(uj,uj',o,o')e l(ul - ul "> t/h A(uj',uj,o',o)dujduj'dodo' (26) 

JO Jo Jo Jo' 
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Considering that A is arbitrary, we obtain the final equation for the evolution of the functional p as: 



Pit) 



p p p OO pQO p p 

/ / p(w,o,o')(w,o,o'\dujdodo' + / / / / p{^,J,o,d)e i ^-^ t l n {uj,J,o,d\dwd^i'dodd (27) 

Jo Jo' Jo Jo Jo Jo' 



where we will call the first term "invariant part" and the second term "fluctuating part" of p(t). 

If we now consider the states p such that the product p(ui, u)' , o, d) A(u>, u>', o, d) is integrable, the Riemann-Lebesgue 
theorem [16] can be applied to eq.(26) to conclude that: 

lim [p{t)\A) = [ [ [ p(u>,o,d)A(u, 0,0') du> dodo' = (p*\A) for any A E A (28) 

Jo Jo Jo' 

where the functional p* is precisely the singular component ps of p (see eq.(13)): 

pOC p p 

p* = / / p(u>,o,o')(iJ,o,o'\cLjdodo' (29) 

Jo Jo Jo' 

The physical meaning of this process can be understood when we consider that the mean value of the observable A 
in the state p can be computed as (A) p = (p\A). Therefore, eq.(28) can be rewritten as: 

lim (A) {t) = (A) p , for any A E A (30) 

Of course, this limit does not contradict the fact that the off-diagonal terms of a functional p representing the quantum 
state of a closed system never vanish through the unitary evolution described by the Liouville-von Neumann equation. 
What self-induced decoherence shows is that the mean value (A) p ^ of any observable A E A will evolve in such 
a way that, for t — > 00, it can be computed as if the system were in a state p* where the off-diagonal terms have 
vanished. Formally this is expressed by the fact that, although we strictly obtain the limit (30) (or (28)), the state 
p(t) has only a weak limit 

w - lim p(t) = p* (31) 

t— >oo 

This weak limit means that, even if p(t) always follows a unitary evolution, the system decoheres from an observational 
point of view, that is, from the viewpoint given by the observable A, for any A E A. 



B. Decoherence in the remaining variables 

As we have seen, for t — > 00 the system decoheres in energy since p* turns out to be diagonal in u>. However, we 
would like to obtain a state diagonal in all the variables. To the extent that we have taken the limit t — > 00, it is 
impossible that a new process diagonalizes the o— variables. As we will see, when a convenient basis is chosen, the 
diagonalization of p can be completed. This second stage necessarily depends on the initial condition p at t = 0, since 
p* is a constant of motion. 

Let us consider a unitary operator U that keeps the Hamiltonian invariant but changes the set of observables 
{Oi, 02, • ■ • , On} into the set {Pi, P2, . ■ . , Pn}, where {H 7 Pi, P 2 , . . . , Pn} is also a CSCO. The simplest form of U is 
given by: 

U= / / U(uj,p, o) \lo,p){lo, o\ duj dp do (32) 
Jo Jp Jo 

The action of U on the ket |w, 0} defines the action of U on any ket, since the kets \ui, o) belong to a generalized basis. 
This action can be easily computed as: 

\w,p) := U\cj,o) = U(u,p,o)\uj,o)dp (33) 
Jp 

The unitarity of U implies that UU^ 1 = I. From here, we obtain: 

U(u,p, o) U*{co,p', o) do = 5(p - p') (34) 



G 



In the new representation, the operator A takes the form: 

A = f (( A(uj,p,p')\uj,p)(uj,p'\dujdpdp'. (35) 

JO Jp Jp 1 

Introducing eq.(33) into eq.(35) gives: 

A= [ ( ( ( ( U{oj,p,o) A(u,p,p')U*{uj,p',o')\u,o)(u,o'\dujdpdp' dodo' (36) 

J J p J p' J o J o' 

Therefore, the coordinates of A in the old basis are: 

A(u>,o,o')= [ ( U(u,p,o) A(cj,p,p')U*(uj,p',o')dpdp' (37) 
J p J p r 

Since U is a unitary operator, eq.(37) is invertible: 

A(u),p,p') = / / U(u),p',o') A(w,o, o')U*(cj,p, o)dodo' 

J O J o' 



(38) 



and, finally, by duality one finds that: 

p(u>,p,p') — / / U(uj 7 p' ,o') p(u>,o, o') U* (oj,p 7 o) dodo' (39) 

J o J o' 

If p = / °° J J o , p(u, o, o') \u, o)(u, o' \ dui dodo' is a state, it must be positively defined and self adjoint. This implies 
that p(ui,o 7 o') = p*(ui,o\o) [13]. Then, we can choose U(oj,p 7 o) such that: 

p(u,p,p') = p(u,p)5(p-p') (40) 

and this completes the diagonalization. In the new basis, eqs.(28) and (29) of the previous subsection become: 

lim(p(t)\A) = [ [ p(u>,p)A(uj,p)dujdp = (p*\A) for any A e A (41) 
'- >0 ° Jp Jo 

P* = p(u,p)(u,p\dLddp (42) 

Jp Jo 

where now p» is completely diagonal in uj and p. The generalized basis given by {\ui,p,p'), \uj,u>' ,p,p')} is the 
preferred basis (also called "pointer basis" 7 ), as presented in [4] and extensively discussed in [5]. On the other hand, 
the decoherence time to can be computed 8 as in paper [15]. 

IV. THE WIGNER TRANSFORMATION 

A. Characterization of the Wigner transformation 

The Weyl transformation maps functions or generalized functions on phase space into operators [1]. Thus, the 
Wigner transformation maps operators into functions on the phase space [17], [18]. If A is an operator, we will 
denote the function corresponding to A via the Wigner transformation by symbA or A(<f>), where <\> = (q, p) = 
(qi, </2, • ■ • , Qn+i,Pi,P2, ■ ■ ■ ,Pn+i) is a point in flat phase space. The function A((p) = symbA is called Wigner symbol 
or Wigner function of the operator A. As the Weyl transformation is a one to one mapping, the image of the algebras 



7 We prefer to use the term "preferred basis" instead of "pointer basis" since it arises not only in measurements contexts where 
the pointer of a measuring device is involved. 
8 Of course, in the case where to — > oo, the system does not decohere. 
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A, As and Ar by the Wigner transformation are non-commutative algebras of functions denoted by C, Cs and Cr 
respectively, where C — Cs®Cr. Since the image of A by the Wigner transformation is the function symbA, it seems 
quite natural to use the notation symb to denote the Wigner transformation itself, so that: 

symb : A i — > C , symb : Ar i — ► £r , symb : As 1 — ► C-s (43) 

Let us define the mapping Ar i — ► £_r for regular observables as usual [18]. First, consider the phase space (in this 
case 9 R 2 ( Ar+1 ^) and endow it with the symplectic form: 

Then, let / be an operator such that symbf = f{4>). This transformation is defined by the usual Wigner recipe as: 

symbf = f{4>) := J d 2JV +Vexp Qv'W 6 ) Tr (t(V0/) (45) 
where ip a and fc denote the a-th and the 6-th components of the points ip and <f> on phase space, respectively. Here: 



where: 



Tty) = cxp ( -r^ h ) (46) 



d d 

<t>=(Qi,-- -,Qn+i, -ih-^— i • ■ • > ) (47) 



and i = 1, 2, . . . , N + 1, is the i-th component of the position operator on the Hilbert space L 2 (R N+1 ). 

The non-commutative product that corresponds to the product of operators in £ (or Cs and Cr) is called star 
product, and it is given by: 

symb(fg) = symbf * symbg = (/ * g)((f>) (48) 

where f(</>) and (/(</>) are the Wigner symbols of the operators / and g, respectively. It can be proven ( [17], eq.(2.59); 
for more general expansions, see [19]) that: 

- /Wexp (yK^ ab ^ 6 ) flW =5^)exp ^-^K^ ab ^ 6 ) (49) 

The Moyal bracket is the Wigner symbol corresponding to the commutator in £: 

{/, g}mb = j^(f*9-9*f) = symb (j^if, g}^j (50) 

Then, if we expand the last two equations in power series of H, we obtain [1]: 

(/ * g)W = f(<f>)g(<f>) + £ h r P r {f{4>)g{<t>)) (51) 

r=l 



r=l 

where the P r are the coefficients obtained by means of eq.(49) and pb means Poisson bracket. This suggests that, 
in the limit h — > 0, the star product should become the ordinary product and the Moyal bracket should become 



9 The fact that the dimension of the phase space is 2(N + 1), where N + 1 is the number of observables of the CSCO, amounts 
to the integrability of the classical system resulting from the Wigner transformation. Non-integrable cases will be considered 
elsewhere. 
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the Poisson bracket. In fact, this is the case in many circumstances although not in all, because in some cases the 
coefficients P r may contain factors of the form h" 1 , making the limit h — > singular (see [1], [3]); in those cases, 
the problem 3 mentioned in the Introduction arises. From eq.(49) we see that factors ft -1 can only come from the 
symbols /(</>) or g(0); then, if these functions do not depend on ft -1 , the limit h — > is regular and can be considered 
as the proper macroscopic limit. 

Finally, let us observe that if / commute with g, eqs.(49) and (51) become: 

(/ * 9){4>) = f{4>) cos (-**d „w°%) g(4>) (53) 
and, hence, in the simple cases with no factors h" 1 we obtain: 

(/*ff)M = /(#K#+o(& 2 ) (54) 



B. The Wigner transformation of observables and states 



In the previous subsection we have considered the Wigner transformation for regular observables, as it is usually 
defined. But the Wigner transformation has not been defined when singular distributions are involved; therefore, the 
transformation must be defined in this case. 

Let us go back to the algebra As of the observables that commute with the CSCO {H, Pi, ... , Pjv}, also denoted 
by {H, P} for simplicity. An element of this algebra is given by: 

As = I / A(u,p) \u,p) ckodp (55) 
Jo Jp 

where A(w,p) is a regular function on its variables. The functional calculus gives: 

A S = A(H,P) (56) 

where P := (Pi, P2, ■ . ■ , P/v). Since we are assuming that the observables H and Pi are Weyl observables, they have 
Wigner functions H{<j>) and Pi{<p), i = 1, 2, . . . , N, respectively, not depending on ft -1 . As a consequence, if / and g 
are any powers of H or Pj, cq.(54) holds. Therefore, if we also assume that A{u,p) is analytic on its variables, we 
have that: 

symbAs = A s (<j>) = A{H {</>), P(tf>)) + 0(h 2 ) (57) 

This means that the problem 3 mentioned in the Introduction does not arise in the singular algebra As when we work 
with Weyl operators. As a consequence, the Wigner symbol of any observable A G As, when h — ► 0, is A{H{<j>), P(4>))- 
Now let us study the Wigner transformation for states of A*. There are two cases that we have to consider. The 
first case is given by the states in A* that can be written as regular density operators. These states are characterized 
by p(u),w,p,p') — p(uj,p,p') (see [7], [8], [10], [12], [13]). Regular density operators have well defined Wigner functions 
[18]: we only must add a (2itK) n+1 factor to eq.(45) in order to obtain the usual normalization for the Wigner function. 
The second case includes any other possibility, i.e., density operators for which either p(ui,w,p,p') ^ p(u),p,p') or 
p(u,u,p,p') or p(w,p,p') are defined only in a distributional sense and not as regular functions. The question arises 
about whether the operators of this second type do or do not have a well defined Wigner function. The answer is 
given by the duality formula (18). If p G ^4* and A G A, we can define symbp = p(<f>) in such a way that it satisfies: 

(symbp\symbA) : = (p\A) = / / / A(u>,p,p') p(uj,p,p') duo dpdp' + 

J a Jp Jp' 

I I jj A(lo,lo' ,p,p') p(u),u>' ,p,p') duj duj' dpdp' (58) 
Jo Jo Jp Jp' 

As the integrals in (58) are well defined, for any p G A*, symbp is also well defined and belongs to the dual space C* 
of the algebra C. Let us recall the decomposition C* = C* s © C* R and the fact that the operation symb is a bijection: 
symbA s 1 — > £%■ Definition (58) will allow us to obtain symbps in the next section. 
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V. THE CLASSICAL LIMIT 



We have seen that, as the result of decoherence, the regular part pn of p vanishes and only the singular part ps = p* 
remains (see eq.(42)): 

P*=Ps= I I p(^,p)(uJ,p\dL;dp (59) 
Jo Jp 

Then, the problem is to find the classical distribution p c {<P) resulting from applying the macroscopic limit h — > to 
the Wigner transformation of p* = ps- 

pOC p pOO p 

Pc{4>) = nm symbps = lim / / p(u,p) symb(uo,p\duj dp = / / p(u>,p) lim. [symb(u> , p\] duj dp (60) 
" h ^Ja Jp Jo Jp 

The first step consists in obtaining the limit of symb\oj,p) for h — > 0. From eqs.(55) and (57) we know that: 

lim [symbAs] = A(H(<f>), P{<j>)) = lim / / A(u,p) symb\u,p) duj dp = / / A(uj,p) [lim symb\uj,pj\ duj dp (61) 
' h ^ a Jo J p ' Jo J P 

But A(H((j)), P(4>)) can also be written as: 

A{H{<j>), P(<j>)) = f f A{u,p) 5(H(<f>) - uj) 5(P(cj>) - p) dpduj (62) 
Jo Jp 



By comparing cq.(61) and eq.(62) we can conclude that 10 : 

lim symb \u),p) = 8(H(<fi) — ui) 5(P{<f>) — p) (63) 

h—>0 

The second step begins by remembering that, from eq.(58) and eq.(ll), symb(uj 7 p\ must satisfy: 

(symb(ui,p\ \symb\w',p')) = {u,p\w',p') = 5(u) - u')S(p-p') (64) 
Since the r.h.s. of the last equation does not depend on ft, the limit for h — > results: 

lim(symb( u>,p\ \symb \ui',p')) = (lim [symb{ u),p\] | lim [symb \u>',p')]) = 5(lu — u>') 5(p — p') (65) 

We will use psio P (<f>) to denote the limit for h — > of symb(u,p\: 

Ps^p{4>) ■= lim symb(u),p\ (66) 
n— >0 

Therefore, replacing eq.(63) and eq.(66) into eq.(65) we obtain: 

(ps^pW | 5{H{4>) - J) <5(P(0) - p')) = 5{u - J) S(p - p>) (67) 

The final step consists in obtaining psuj P {<j>)- Since we have assumed that the number of operators in 
the CSCO {H, P\, . . . , Pjv} coincides with the number of degrees of freedom of the system under consider- 
ation and, as a consecuence, the phase space has dimension 2{N + 1) (see footnote 9), then there exists 
a canonical transformation that carries the position-momentum variables <fi = (<L P) m to the variables ip = 
(T(0),ai(0),...,a JV (0),fl"(^),Pi(^),...,Pj V (0)) ) where H(<f>) = symbH, P(0) = symbP h i = 1,2,. ..,N, r(0) is 
the conjugate variable of H(<f>), and the a,(0) are the conjugate variables of the Pi (</>). On the other hand, from 
subsection III. A we know that all the singular operators As £ As and all the singular functionals ps £ A* s are time 
invariant. Since the function symb does not introduce time variables, all the As{4>) = symb As £ C s and all the 



10 In the case of discrete spectrum we would have A(u),p) = S(lo — Lo')8 pp , instead of A(u>,p) = S(lu — to 1 ) S(p — p'). Then, we 
would obtain symb\u>' ,p) = 8(H(4>) — u')Sp^ p i, where Sp^ p , is a N Kronecker 5. 
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ps(<P) — symbps £ C* s are also time invariant, and the same holds for their limits for h — > 0. Thus, if A{<j>) £ £5, 
p{4>) £ jC* s , A{4>) = f(H(4>),P(<t>)), and p{4>) = g(H (<$>), P{4>)), then (p(<j))\A(<j))) can be expressed as: 

( P {4>)\A{4>)) = [ f(H(<t>),P(<j>)) g(H(<j>),P(cP)) # 2 ( JV + 1 ) = fdr [ da [ f f(H,P) g(H,P)dHdP (68) 
J02(jv+i) j T j a JhJp 

We can call: 

C{H, P)= [ dr da (69) 

JM(H,P) 

where C(H,P) is the volume of the region A4(H,P) of the configuration manifold defined by the conditions H = 
const, and P — const. Then, eq.(68) becomes: 

(p(<f>)\A(<t>))=C(H,P) f [ f(H,P)g(H,P)dHdP (70) 
Jh Jp 

In the simplest case of bounded integrable systems described by action-angle variables, C (H, P) is a constant equal 
to (2ir) N+1 . Anyway, C(H,P) is always a constant that we will ignore from now on in order to simplify notation. If 
we now apply the result expressed by eq.(70) to eq.(67), we obtain: 

/ / Psu P {H,P)5(H -uj')S(P - P ')dHdP = 8{uj-J)5{p-p') (71) 
Jh Jp 

This means that 11 : 

psu P {H,P)=6{H-uj)5(P-p) (72) 

If we now go back to the variables </>: 

p s „p(<j>) - lim symb(Lu,p\ - 8{H{cj>) - u) 5(P(<f>) - p) (73) 
h—>o 

Finally, we can obtain the classical distribution p c by replacing the just obtained result (73) into eq.(60): 

Pc{</>)= [ f p{u,p)psup{4>)du)dp= [ [ p(uj,p)S(H((f>)-uj)S(P{<p)-p)dujdp (74) 
Jo Jp Jo Jp 

As we can see, p c {4>) is a constant of motion, as it was expected due to the results obtained in Subsection III. A. 
Eq.(74) has a clear physical meaning: p c {4>) is a sum of densities infinitely strongly peaked on the classical trajectories 
defined by the constants of motion H{<j>) = u and P{4>) = p and averaged by the density function p(u),p) which is 
properly normalized according to eq.(16). As a consequence, p c {4>) can be conceived as sum of classical trajectories 
weighted by their corresponding probabilities. This leads to the expected result: classical motion takes place along a 
classical trajectory, and the probability of each possible trajectory is given by the initial condition p at t = 0. 

It is also interesting to consider the case where the initial condition p at t — is such that the factor 
p(u),u>' ,p,p') A(u),u>' ,p,p') of eq.(26) is strongly peaked around u> — oj' = E. In this case, the evolution factor 
e i{tt}-tt) )t/h can approximated by e lEt / h . This shows two facts. First, there is an interplay between the character- 
istic energy E and the decoherence time: the decoherence time becomes shorter as the energy is higher. Second, the 
limit E — > 00 plays the same mathematical role as t — > 00, and represents the well known "high energy limit" [3]: for 



11 In the discrete case, eq.(71) reads: 

/ / Ps„ P {<t>)5(H - w')S^ p ,dHdP N = S(u-u')S. 
Jh J p 

Then: 

Asu,V W = 8{H{4>) - cu') 5 N {P(4) - p') 
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high energies many systems behave in an almost classical way, e.g., high energy orbits of atoms can be approximated 
by classical trajectories 12 . 

In this section we have found the classical limit by applying the macroscopic limit to the Wigner transformation of 
the state resulting from decoherence. Nevertheless, it is also possible to translate the quantum evolution equation 
into the classical language via the Wigner transformation in order to obtain the complete process of decoherence in 
classical terms. This strategy leads to the same result as the one obtained in the present section (see Appendix A). 

Finally, let us recall the three problems of the traditional way of conceiving the classical limit as they were presented 
in the Introduction, and let us consider how and under what conditions they can be overcome from the present 
approach: 

1. Although in general the Wigner function p(<j>) of a state p is not non- negatively defined, we can guarantee the 
non-negativeness of p c {4>)- I n fact, p(u>,p) is non-negatively defined due to its origin, since it represents the 
diagonal components of the original state p (it can be also proved that, in general, if p = p(u>) = p(u>,u>), then 
the limit for H — > of p(<p) is non-negatively defined a.e.; see Appendix B). 

2. Although only Hamiltonians of degree < 2 in p and q yield to Hamiltonian fluxes that maintain the deformation 
invariant, this is not a problem from the present perspective since, after decoherence, only the singular algebra 
remains, and in this algebra both states and observables are constants of motion. 

3. Although in some cases factors of the form hr 1 may appear in the Wigner state function making the limit h — > 
singular, we have shown that this possibility is blocked when the observables of the CSCO are Weyl observables. 
The requirement of working with a CSCO consisting of Weyl observables is very weak since it does not impose 
artificial constraints on the state p. Furthermore, to apply the Weyl transformation to classical observables is 
the usual strategy in practice for obtaining the corresponding quantum observables. 

VI. THE PHYSICAL MEANING OF THE CLASSICAL LIMIT 

As we have seen, the classical limit of quantum mechanics involves two elements: 

1. Decoherence: According to the self-induced approach, decoherence is a physical process that depends on the 
own dynamics of a closed quantum system governed by a Hamiltonian with continuous spectrum. As the result 
of decoherence, in the infinite time limit the mean value of any relevant observable can be computed as if the 
system were in the diagonal state /?*. In other words, decoherence transforms standard (non-Boolean) quantum 
mechanics into a Boolean quantum mechanics restricted to states that are diagonal in the basis defined by the 
CSCO {H,P}. 

2. Macroscopicity: For ft-»0, the Wigner transformation of the diagonal state p* turns out to be p c ((f>), and it is 
resolved into an ensemble of classical trajectories on phase space weighted by their corresponding probabilities. 
This means that, in the macroscopic limit, the Wigner transformation maps the Boolean description resulting 
from decoherence onto classical statistical mechanics. 

This shows that, strictly speaking, the classical limit of quantum mechanics is not classical mechanics but classical 
statistical mechanics. This point deserves some further remarks. 

In the classical distribution p c (4>) resulting from the classical limit, the ensemble of trajectories is weighted by 
the non-negative function p(u>,p): it is precisely the fact that p{iv,p) is non-negatively defined what permits it to 
be interpreted as a probability function. But the formal agreement between p c (4>) and a density distribution in 
standard classical statistical mechanics does not mean that both have the same physical meaning. In fact, in classical 
statistical mechanics probabilities are conceived as a sort of measure of our ignorance about the real deterministic 
classical trajectory. On the contrary, since the p(uj,p) are the diagonal components of the original quantum state 
p, they represent quantum probabilities which, as many no-go theorems show, are irreducible. Of course, this fact 
does not mean that a particular classical trajectory cannot be chosen. Let us suppose that we prepare the system at 
t = in an initial condition p such that its singular part ps is an eigenstate (u),p\. In this case, as a consequence of 



It is worth noting, however, that while t is a perfectly well defined variable, E is just a characteristic energy that may be 
not well defined in some cases. For this reason, the mathematically precise strategy is to find the limit for t — ► oo as we have 
done in the present work. 
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decoherence and macroscopicity, we will obtain the particular trajectory defined by the constants of motion H{<j>) = lo 
and P{4>) — p with certainty. This shows that, although p(u),p) in the classical distribution p c {4>) represents quantum 
irreducible probabilities, a particular classical trajectory can always be selected by means of the proper preparation 
of the quantum initial conditions. 

Finally, it is worth stressing the emergent nature of classicality as explained by the present approach. As we have 
seen, the off-diagonal terms of the quantum state p(t) never vanish through the unitary quantum evolution. Strictly 
speaking, what self-induced approach shows is that, in the infinite time limit, for any A G A , {A) p ( t ) can be computed 
as if the system were in the diagonal state p*. In fact, the limit for t — > oo of (A) p u\ could also be computed in the 
Hciscnberg picture as the limit for t — > oo of (A(t)) p ; in this case we would obtain a diagonal operator ^4*. This fact 
clearly shows that the fundamental magnitude in the explanation of decoherence is {A) p ( t ) = (A(t)) p and not nor 
A*. In other words, decoherence should be conceived as a coarse-grained process that describes the evolution of the 
state p(t) from the observational viewpoint given by the observable A 13, . As a consequence, classicality is an emergent 
property that arises in a coarse-grained level of description . The classical limit shows that, from the point of view 
given by A, as the result of decoherence and macroscopicity the quantum system behaves as if it were a classical 
statistical system. This means that our measurements of the mean value of any relevant observable A on the quantum 
system will give the same results as those we would obtain on a classical statistical system described as an ensemble 
of classical trajectories weighted by their corresponding probabilities. The distinction between the fundamental and 
the coarse-grained levels of description permits us to understand how the Boolean and deterministic classical world 
objectively arises from an underlying non-Boolean and indeterministic quantum level. 

VII. CONCLUSION 

Einstein was right when he considered the idea that the limit h — ► is the right classical limit as an oversimplification. 
On the basis of the assumption that the problem of the classical limit of quantum mechanics amounts to the question of 
how the classical world arises from an underlying quantum reality, our account of the problem involves two elements. 
The first one is self-induced decoherence, conceived as a process that depends on the own dynamics of a closed 
quantum system governed by a Hamiltonian with continuous spectrum; the study of decoherence was addressed by 
means of formal tools derived from the van Hove formalism. The second element is macroscopicity represented by the 
limit h — > 0; we have shown that, when the macroscopic limit is applied to the Wigner transformation of the diagonal 
state resulting from decoherence, the description of the quantum system becomes equivalent to the description of 
an ensemble of classical trajectories on phase space weighted by their corresponding probabilities. Furthermore, this 
approach to the classical limit explains under what conditions the problems arising from the traditional approach can 
be avoided. Finally, when these formal results are considered in the light of a generalized concept of coarse-graining, 
decoherence turns out to be a coarse-grained process that, in the infinite time limit, leads to classicality when the 
system is macroscopic enough. Since there is no subjective element involved in this process, from our approach 
classicality is a property that objectively emerges from the underlying quantum world. 
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13 '{A) p (t) = (p\A) can be thought as a generalized "projection" of the state p. In fact, we can define a projector II belonging 
to A® A* as II = ApA, where pA satisfies (pa\A) — 1 (this condition guarantees that II 2 =11). In this case, p re i — (p\A)pA, 
where p re i is the projected part of p relevant for decoherence. Since coarse-graining amounts to a projection (see [20]), {A) p ^ 
can be conceived as a coarse-grained magnitude. 

14 An interesting discussion about emergence and reductionistic relations between the various levels of the quantum mechanical 
descriptions can be found in papers [21], [22]. 
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APPENDIX A: THE QUANTUM EVOLUTION IN CLASSICAL TERMS 



The quantum evolution of the system can be translated into the classical language via the Wigner transformation. 
The phase space analogue of the Liouville-von Neumann equation for ft — > 0, eq.(52), is: 

d tP {t) = {H,p} mb (Al) 

Let us now compute p((f>, t) = symbp(t), which is equal, up to the order ft 2 (hence, it is in fact the limit when ft — > 0), 
to: 

/■ GO pOO p p 

p(<j>,t)=p(H(cj>),P(ct>)) + / p^^^e^-^^symb^^^p^dwdw'dpdp' (A2) 

JO JO JpJp' 

Here we have used the Liouvillc equation for the regular part, and we have kept the singular part unchanged since we 
know that it is time invariant. As in the quantum case, we can call the first term the "invariant part" and the second 
term the "fluctuating part" of p((j>,t) [10], [15]. On this basis we can also compute: 

(p(<j>,t)\A(<t>))= [ f P (H(<j>),P(cP))A(H(<j>),P(cf>))dHdP + 

J H J P 

pOC pOC p p 

/ / / / piu^^p,^)^^-"'^^ A{u,us\p,p')dujduj' dpdp' (A3) 

Jo JO Jp Jp> 

where, again, we call the first term "invariant part" and the second term "fluctuating part" . Now, if the product 
p(u) , u)' , p, p')A{lo' , u , p' , p) is integrable, we can use the Riemann-Lebesgue theorem to conclude that: 

\im (p(<t>,t)\A(<j>))= [ [ p(H(4>),P{4>)) A(H(4>),P{4>)) dH dP = {p*{<f>)\A{4>)) (A4) 
JhJp 

for any A(<f>) G L and for any p(<p) G C* with the right properties. Therefore: 

P* (</>)= / / P {H{4,),P{4>))dHdP = f°° f p(u,p)5(H(cP)-w)5(P( ( l>)-p)ckjdp (A5) 

JH JP Jo Jp 

This last equation is precisely the Wigner transformation of for h — > 0, as obtained in eq.(74). In p*((f)) the 
non-diagonal terms have disappeared and only the diagonal (singular) terms remain. Thus, we have found the weak 
limit: 



w - lim p{<p,t) = p*(4>) (A6) 

t— »oo 



which express the result of decoherence. 



APPENDIX B: POSITIVITY OF THE WIGNER FUNCTION OF p(u) 

Here we will prove that the Wigner function of a state represented by a density operator of the form p(u) = p(u, to) 
is positively defined a.e. in the limit fi->0. This proof is a reformulation of an argument due to Narcovich [23], [24]. 
Let us call pn(q,p) the Wigner function for the density operator p in terms of h. We will call: 

Jim Ph(q,p) = G(q,p) (Bl) 
We will prove that G(q i p) is non-negative (a.e.): 

G(q,p) > (B2) 

Let us call: 

•-(;)■ 2= 0- '=('3 <B3) 

14 



and 

a{a, z) = aJz = (q,p) ( ] f ? ) = gp' - p?' (B4) 

Now we consider the inverse symplectic Fourier transformation of pn(q,p) = pnifi) given by: 

h(z) = J pnWeW^dz (B5) 

where dz = dqdp. This transformation (i.e. the function fn{z)) is of ^-positive type (see [23], [24]), which means 
that: 

m 

]T fh(aj - o fc ) e 4(/l/2) *<°*'°'> X*X k > (B6) 

where: 

1. - a 1; a 2 , . . . , a m are arbitrary points of the phase space: 

a k = I q j (B7) 

2. - Ai, A 2 , . . . , \ m are arbitrary numbers. 

3. - m = 1,2,... is an arbitrary finite positive integer. 
Let us now consider: 

9{z) = \vcah(z) (B8) 

h—>o 

and see whether g(z) is positive in the sense of Bochner. If we make the limit K — > in cq.(B5), we obtain: 

9(aj - a k ) \*\ k > (B9) 

j,k=l 

where the notation has been defined in 1, 2, 3 right above. Thus, the function g(z) is positive in the sense of Bochner. 
Therefore, the Fourier transform p(q,p) is also positive (a.e.): 

¥>(«,P) = J g(z)e~* za dz>0 (BIO) 

This property is inherited by the symplectic Fourier transform of g(z), that we have called G(q,p): 

G(q,p) = <p(p,-q) (Bll) 

Then, from eq.(B8) we have: 

lim p h (q,p) = f [lim f H {z)] dz 
hi— >o j ih^o 

= j g{z)e- ll7 ^ z) dz = G(q,p) > (a.e) (B12) 

where pn(q,p) is an arbitrary regular Wigner function, q.e.d. Nevertheless, it is worth stressing that the non- 
negativeness obtained through decoherence is stronger than this result: p c {<t>) is non- negatively defined on the whole 
phase space, and not only a.e. 
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